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Abstract 

We develop the general formalism of string scattering from decaying D-branes in 
bosonic string theory. In worldsheet perturbation theory, amplitudes can be written as a 
sum of correlators in a grand canonical ensemble of unitary random matrix models, with 
time setting the fugacity. An approach employed in the past for computing amplitudes in 
this theory involves an unjustified analytic continuation from special integer momenta. We 
give an alternative formulation which is well-defined for general momenta. We study the 
emission of closed strings from a decaying D-brane with initial conditions perturbed by the 
addition of an open string vertex operator. Using an integral formula due to Selberg, the 
relevant amplitude is expressed in closed form in terms of zeta functions. Perturbing the 
initial state can suppress or enhance the emission of high energy closed strings for extended 
branes, but enhances it for DO-branes. The closed string two point function is expressed 
as a sum of Toeplitz determinants of certain hypergeometric functions. A large N limit 
theorem due to Szego, and its extension due to Borodin and Okounkov, permits us to 
compute approximate results showing that previous naive analytic continuations amount 
to the large N approximation of the full result. We also give a free fermion formulation of 
scattering from decaying D-branes and describe the relation to a grand canonical ensemble 
for a 2d Coulomb gas. 
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1. Introduction 



D-branes, as solitons of open string theory that are localized in space, have given many 
insights into nonperturbative phenomena in string theory such as string duality, and resolve 
many timelike singularities of General Relativity including those inside some black holes. 
In addition, they lead to the holographic description of 11 dimensional asymptotically flat 
space via a Matrix model and of spaces with a negative cosmological constant in terms of 
a dual conformal field theory. Many of these developments, and particularly the last, arose 
from an understanding of D-brane dynamics - namely how closed strings scatter from such 
solitons which are quantized in terms of open string fluctuations. 

It is of great interest to understand similar issues in the time dependent context of 
rapidly expanding universes, particularly in view of the likely occurrence of a Big Bang 
followed by inflation in the early universe and the possible presence of a positive cosmo- 
logical constant A > now. Exploration of the symmetries and structure of universes 
with positive A has suggested that if they have a holographic description the dual might 
be related to a Euclidean CFT living on the early or late time boundaries of such space- 
times (i, i)i Time, in such a picture, emerges holographically via the RG flow of the 
Euclidean field theory dual in analogy with the emergence of the radial direction of AdS 
space from the RG flow of a Lorentzian field theory (0, ||). In order for such a picture 
to be actually realized in string theory one would need some kind of D-brane localized in 
time, called an S-brane in ||, with a decoupling limit relating the closed strings on the 
spacetime background to the open strings quantizing the brane. One might hope that 
the spacetime near such a Euclidean brane would be rapidly expanding by analogy with 
the rapid expansion of the transverse space in the vicinity of a conventional brane. In 
the decoupling limit for standard D-branes it is precisely this rapid transverse expansion 
that gives rise to the exponential increase in the volume of AdS spaces as the boundary 
is approached. Thus one might hope that a Euclidean brane would lead to a spacetime 
exponentially expanding in time in a suitable decoupling limit. 

Sen has proposed |Q] that S-branes are concretely realized in string theory by the 
exact boundary CFTs representing decaying branes. In bosonic string theory the exactly 
marginal boundary interaction 



(where X° is the timelike scalar) describes such a brane ||, 0, The basic results 

that have been obtained so far are: 
1. At vanishing g s the brane decays to "tachyon matter" with energy but vanishing 
pressure [Q. 

5 One heuristic way of motivating this is to note that Lorentzian de Sitter space and Euclidean 
AdS space are solutions to — Xq + Xf + • • • X\ = ±l 2 embedded in flat (d + 1) dimensional space. 
Both hyperboloids are thus different real sections of the same complex manifold. 

6 There are also coshX and sinhX interactions, which describe brane formation and subse- 
quent decay [Q. Note the analogy with the inflationary and global parametrizations of de Sitter 
space, namely ds = —dt + e dx versus ds = —dt + cosh(t) dil . 
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2. At finite coupling there is tremendous production of very heavy non-relativistic closed 
strings 

3. It has been proposed that the resulting coherent state of heavy closed strings is an 
equivalent description of the tachyon matter. 

4. In view of (2) and (3) it is proposed that there is a new open-closed duality hinting 
at a new kind of holography , |10| . 



5. There is a close relation with the dynamics of decaying unstable branes in the c = 1 
matrix model that suggests that the picture (1) - (4) is essentially correct (||11||, jl2 

BE 



To understand the structure of tachyon matter, and to explore the possibility of a 
decoupling limit leading to timelike holography, a central problem is to compute general 
scattering amplitudes of closed and open strings from the decaying brane. Closed strings 
can be used particularly to probe the structure of the tachyon matter final state that exists 
after the brane decays, and open strings, since they will only be present in the spectrum 
at early times, can be used to change the initial brane configuration which is decaying. 

Thus, in Sec. 2 and 3 we develop the general formalism of string scattering from de- 
caying branes in bosonic string theory. In a perturbative approach, amplitudes can be 
written as a sum of correlators in a grand canonical ensemble of unitary matrix models, 
with time setting the fugacity. If we integrate over the zero-mode time coordinate at the 
outset, the amplitudes are connected by analytic continuation to a single unitary matrix 
model, the rank of which is related to the total energy carried by the external vertex op- 
erators. Based on the approach which was employed with success in bulk Liouville theory 
|T5|[ , it was proposed in (0, [|16]) to compute bulk correlation functions by "analytically 
continuing" the vertex operator momenta from discrete integer values where direct com- 
putation is simple. This procedure has long been known to be somewhat questionable 
since an analytic function cannot be defined from just the data at a discrete set of points 
without additional constraints such as knowledge of the behavior at infinity or consistency 
conditions imposed by symmetries. In the case of Liouville theory, suitable consistency 
conditions were found (see the review ||17|| ) and one can try to exploit the resulting tech- 
niques to explore the physics of decaying branes (0, [fTB| , |Ji~9fl , pOj] , pl|). For full branes, 



an interesting prescription to compute n-point closed string disk amplitudes was given in 
terms of viewing the brane as an array of ordinary D-branes in imaginary time ||22|| (for 
generalization to higher genus, see |[23|| ). In this paper, we describe an alternative formula- 
tion of the scattering amplitude calculation which is well-defined in an open subset of the 
complex energy plane and therefore permits reliable analytic continuation. We will study 
in detail the open-closed and closed-closed two point functions, and compare our results 
to previous work. 

In Sec. 4, we study the emission of closed strings from a decaying brane with initial 
conditions perturbed by the addition of an open string vertex operator .0 Using an integral 
formula due to Selberg, the relevant amplitude is expressed in closed form in terms of 
zeta functions. Perturbing the initial state in this way can either suppress or enhance 
the emission of high energy closed strings for extended branes, but always enhances it for 



Another approach to these amplitudes will be described in [21]. 
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DO-branes. This is consistent with the picture that DO-branes decay entirely into closed 
strings, with very heavy closed strings making up the tachyon matter state, while higher 
dimensional branes can decay partly into unstable lower dimensional branes. However, the 
enhancement of emission that we find in some cases will increase the direct closed string 
emission even by the higher dimensional branes. We discuss the consequences for the 
hypothesis that tachyon matter is nothing but a state of very heavy closed strings. In Sec. 5, 
the closed string two point function is expressed as a sum on N of Toeplitz determinants of 
certain NxN matrices of hypergeometric functions. A large-N limit theorem due to Szego, 
and its extension due to Borodin and Okounkov, permit us to compute approximate results 
that show that previously used methods to compute scattering amplitudes in the decaying 
brane amount to a leading large N approximation from our perspective. In Appendix C 
we give a free fermion formulation of scattering from decaying branes by extending old 
techniques of Douglas []24] . In Appendix D we describe the relation of the decaying brane 
correlators to a grand canonical ensemble for a classical Coulomb gas. 



2. Review of closed string tachyon scattering from D-branes 

To establish notation and conventions, in this section we review the standard compu- 
tation of closed string tachyon scattering from a D-brane. This amounts to computing the 
bulk two point function on the disk, with Neumann and Dirichlet boundary conditions in 
the various directions. 

2.1. Kinematics 



We follow the conventions of Polchinski [[25), with r)^ = (— , +, +,+,...), and a' = 1. 
The closed string tachyon vertex operator is 

e lk - x , k 2 = -m 2 = 4. (2.1) 

For scattering from a Dp-brane, we divide the momenta into the p + 1 parallel components 
fe", and the 25 — p transverse components k^. We write the momenta of the two closed 
string tachyons as fci = (k^k^) and k 2 = ( — k "jfeg"). The Mandelstam variables are 
defined as 

s = 2(/cll) 2 , t = k 1 -k 2 . (2.2) 

We can f actor ize the amplitude in the closed or open string channels, and find poles 
at the location of physical string states: 

closed : k 2 = -m 2 = -4(N - 1), 

(2.3) 

open : k 2 = —m 2 = —(N — 1) 
with N = 0, 1, 2, . . .. This implies that poles occur at 

closed : t = -2(N + 1) 

(2.4) 

open : s = -2(N - 1). 
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2.2. Correlators on disk 

We work on the unit disk, \z\ < 1. As is standard, we separate out the zero modes 
from the worldsheet fields X fI (z,z) by writing 

X»(z,z) =x» + x'»(z,z) (2.5) 

with Jd 2 zX >*(z,z) = 0. The zero mode integrals are done at the end of the calculation, 
and enforce momentum conservation in Neumann directions. 
The Neumann and Dirichlet correlators on the disk are! 

-if [ln\z - w\ 2 + zw\ 2 ] N 

(X'»{z,z)X'»{w,w))=l \ (2.6) 

-rf v [ln|z - w\ 2 - In |1 - zw\ 2 ] D. 

For z in the bulk and w on the boundary {w = e lt ): 

(X-{z,z)X"(t)) = l q d (2.7) 

For both points on the boundary: 

-77^ hale" 1 - e it2 \ 2 N 
(Jf>(t 1 )X-"(t 2 )> = < (2.8) 

D. 



For self-contractions in the bulk: 



-rf v - zz\ 2 N 



(X»(z,z)X»(z,z)) = { 2 (2.9) 

-rj^ In \1 - zz\ 2 D. 



For self-contractions on the boundary: 



N 

(X'»(t)X v (t)) = { (2.10) 
D. 



8 In fact, there is a subtlety in defining the Neumann correlator. The Green's function 
which obeys Neumann boundary conditions and vanishes when integrated over the disk is 
— \rf' v [in \z — w\ 2 + In 1 1 — zw\ 2 — zz — ww — c] , where c is a number which depends on the 
choice of worldsheet metric. The extra terms drop out after using spacetime momentum con- 
servation (see section 6.2 of f25|| ), and so don't contribute to scattering from an ordinary D-brane. 
However, such terms do seem to contribute when considering D-branes with nontrivial worldvol- 
ume fields, as is our interest. Nevertheless, following standard practice, we will continue to drop 
these terms by adopting this as our prescription for defining the (naively divergent) X° path 
integral. 
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2.3. Closed string tachyon scattering from static D-brane 

Using the above correlators, the two point function of bulk exponentials is 

( e ifc 1 .*'(z,S) e *fc 2 .X / (ti;,fi>)) = \ z - w \t\\ - zw^s-^l - zz\i- 2 \l - WW\%- 2 . (2.11) 

The S-matrix amplitude is given by fixing one vertex operator at the origin, z = 0, and 
integrating w over the disk. Up to an overall constant and a momentum conserving delta 
function, the amplitude is 

Cdr r r*(l - r 2 )§" 2 = F + *) r (f ~ *) (212) 

Jo r (§ + 1) 

which indeed exhibits poles in accord with (|2.4|) . 

In the remainder of this work we are interested in computing the scattering amplitude 
from a decaying D-brane. Physically, as a function of x , we expect to find an amplitude 
which interpolates between that of a D-brane and that of a collection of closed strings 
(or tachyon matter) into which the D-brane decays. One signature of this is that the 
open string poles should be absent at late times. Besides this, it is difficult to exhibit 
any "smoking gun" signatures of the presence/absence of the D-brane. For instance, it is 
known that in the high energy and fixed angle regime, both types of amplitudes behave as 
e~ E indicating softness at short distance. 



3. Scattering from the rolling tachyon 

Now we consider the tachyon profile T(X°) = Xe x ° corresponding to the boundary 
interaction 

e S bndy =e ~\e*° fdte x ° ( 3 _-g 



where we have separated out the zero mode as in (|2.5|) . Consider the scattering of £ bulk 
tachyons. The X° part of the computation involves 



A e = DX°e- s H 



e -S 77 e iwaX°(z a ,z a ) 



e (3-2) 



= [dx e™°^(e- Xe *° J dteX '° f[e* 

J a=l 



iio a X (z a ,z a ) 



The full amplitude also contains terms involving the spacelike X 1 that are the same as for 
a standard D-brane, as well as integrals over the locations of the vertex operators z a . As 
we will describe below the bulk vertex operators can be moved to the boundary (\z\ = 1), 
to compute amplitudes with arbitrary numbers of open string tachyon vertex operators. 
We can always choose a gauge in which our vertex operators (assuming they carry nonzero 



energy) contain no timelike oscillators [26]. Thus the interesting part of any correlation 
function (e.g., for closed string scattering from a brane with a perturbed initial state) 
involves interactions with the boundary tachyon that are summarized by (|3.2|). 
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3.1. Perturbative approach and matrix integral formulation 



One approach is to expand (|3.2| ) in powers of the boundary interaction; i.e. as a power 
series in A. The magnitude of A can be changed by shifting x , so truncating the power 
series is only sensible when considering scattering processes which are dominated at early 
times. In general, we must keep and sum the entire series. For the cases we consider, the 
sum will have a finite radius of convergence. 

Expanding, we obtain 



/oo 
dx°e ix °^a^ £ 

7V=0 



(-27rAe a 



N\ 



N 



dti 



TT— ( 



,x °(ti) 



,x'°(t N ) 



1=1 



n 

0=1 



,w a X °(z a ,z a )^ 



(3.3) 

By separating out the x° integral in this way we are isolating the contribution to the 
total scattering amplitude from the partially decayed state of brane that is present at any 
particular time. The late time contribution from x° — > 00 should isolate the effects of the 
tachyon matter final state to which the brane is supposed to decay. 
To calculate the fixed x° contributions we now need to evaluate 



a 



(TV) 



N 



2tt 



,x °(t N ) 



e" . . .e" v " J,y Y\_ e 

j=l a=l 



iuj a X °(z a ,z a ) 



) ■ 



(3.4) 



The Wick contractions are straightforwardly evaluated using the Green's functions in (|2 
|), yielding 



a 



(AO 



.a<b 



Z a Zb\ 



ah 



N 

n 



dti 
2tt 



n 



it it • 1 2 

q 1 e 



z n .e 



-iti 1 2itu a 



(3.5) 



An elegant way of rewriting ( |3.5| ) is in terms of matrix integrals []T6| . In particular, the 

measure Yli=i ^ Tii<j \ elti —e lti | 2 is nothing but the measure for U(N) matrices expressed 
in the eigenvalue basis, the product of exponentials being the Vandermonde determinant. 
This leads to the identification 



a[ N) = N\ 



I I \z a ~ Z h \ 

.a<b 



-U) a U)b 



/ ^TT|det(l-^[/)| 2ia 

JU(N) 



(3.6) 

where is U is a unitary NxN matrix, and we have normalized the measure to Ju^\dU = 1. 
We can now write 



.4, 



n 

.a<b 



Z a ~ Z h \ 



-UJ a UJb 



ab 



Jdx°e lx °^a^F(z 1: 



(3.7) 
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where we have isolated the relevant summation by defining 



oo „ 

F(z u wi, . . . , zi, ur, (j)=J2 e~ Nf * / dU TT | det(l 



Z a U)\ 



2iuj a 



(3.J 



and 



-x 



ln(-27rA) 



(3.9) 



Here \i plays the role of a complex chemical potential, controlling the weights of the 
various U(N) matrix integrals. As we will discuss in Appendix D there is indeed a precise 
correspondence to the grand canonical ensemble of a statistical mechanical system. Note 
that the late time behavior corresponds to Re(fx) — > — oo. 

The essential problem is to compute the matrix integrals appearing in ( |3.8|) , perform 
the sum over N, and evaluate the x° integral. Finally, to compute a string theory S-matrix 
element, we must also integrate over the positions of the vertex operators. 

Having obtained ( |3.7| ) for general bulk tachyon amplitudes, it is simple to include 
boundary tachyons. Simply take any number of vertex operators to the boundary, z a — > 1, 

and remove the corresponding factors of |1 — z a z a 
contractions according to ( |2.10|) . 



1 2 

i W. 



2 Ua , since there are no boundary self- 



3.2. Integrated approach: performing the zero mode integral first 

An alternative approach is based on a strategy which was employed with success 
in bulk Liouville theory [15]. Instead of leaving the x° integration until the end, we can 
perform it at the outset. This approach requires analytic continuation in the momenta and 
can miss non-analytic pieces that are known to be present in the amplitude. Nevertheless, 
useful information can be obtained. 

So let us return to ( |3.2| ) and perform the x° integral using 



dx° e~ aeX e luJX 



(3.10) 



to obtain 



X u (t) 



£. — J a 



lb 

a=l 



,iuj a X °(z a ,z a ) 



Using the identity T(l — iz)T{iz) = —in/ sinh-Trz we can write this as 

(27rA)" l E Q ^ 



A 



mt 



-Z7T 



Bp 



sinh(7r^ a w a ) 
1 



r(i-<Ea«-) 



dt_ x \t) 
2tt 



n 

a=l 



ito a X °(z a ,z a ) 



(3.11) 



(3.12) 
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As we will see, the overall factor in A 1 ^ is of exactly the form expected for the closed string 
one-point function ||. The expression for Bg contains complex powers of the fields and 
therefore must be defined by appropriate analytic continuation. In particular, to apply 
standard techniques, we need —i J2 a u a to be a positive integer. We will now follow the 
procedure of defining Bi by continuing the external momenta to imaginary integer values 
and then "continuing" back. This procedure requires some care, since in general the values 
of a function at a discrete set of arguments does not determine its behavior in the entire 
complex plane even with an assumption of analyticity. However, we will find situations in 
which the expression for Be, evaluated formally for imaginary integer momenta, is actually 
defined on some open subset of the complex to plane, thus permitting reliable analytic 
continuation. 

Proceeding in this fashion, we perform the continuation 



lUJr, 



-Tin 



to no n- negative integers n a , finding 



Bi = m 



N 

n 

3=1 



2tt 



n 

a=l 



1 a X °(z a ,z a ) 



(3.13) 



(3.14) 



with N = n a . The correlator is the same as in ( |3.4|) after the substitution ( |3.13| ). The 
result is therefore 



B, 



n 

.a<b 



Zb 



n a n b 



ab 



[ dUT\\det(l -z a U)\~ 2na . (3.15) 

JU(N) 



Comparing with ( |3.7| ) an essential difference is that there is no longer a sum over N which 
needs to be evaluated, and the exponents in ( |3.15|) are restricted to special integer values 
(which we'll eventually have to continue to generic values). 

At these imaginary integer momenta, and for sufficiently large N, the matrix integral 
( |3.15|) is iV independent, and can be efficiently evaluated group theoretically |16 



or using 



SU(2) current algebra 0. As shown in |16j, for iV > J2 a n a (which applies to our case, 
since N = Y. a n a), 



/"dC/J]|det(l-^C/)|- 2 ^ = nil 

a ab 



Z a Zb 



-n a n b 



(3.16) 



This then yields 

Be = J] I 2 - - z b\ nan " II I 1 ~ *aZb\~* nanb - (3.17) 

a<b ab 

The final step is to "undo" the analytic continuation, by replacing n a = —iu) a in 
( [3. 171) , which, as we have stressed, cannot be justified without additional information. 
Nevertheless, if we optimistically apply this procedure we are led to: 



/lint 



—17T 



sinh(7rj] a w a ) 



HlZa-Zbl-^llll-ZaZbll 



U! a U! b 



(3.18) 



a<b 
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In Sec. 4 we will see that naively using Q3.16 ) and then analytically continuing to real 
momenta leads to unphysical results for bulk-boundary amplitudes and we will describe a 
better procedure for computing Bi. In Sec. 5 we discuss different procedures for computing 
bulk-bulk amplitudes. 



3. 3. Comparison of one-point functions 

It is helpful to compare the results for the tachyon one-point function (£ = 1) obtained 
in the two approaches.^ By a conformal transformation we can take the bulk vertex 
operator to be at the origin of the disk, z = 0. In the perturbative approach we find from 
([Df) that F = 1/(1 + 27rAe x °), and therefore 

Ax = / dx° = -Z7r (27r ^~" . (3.19) 

/ l + 2n\e x ° sinhTrw v ' 



Comparing with (|3.18| ) we then find an agreement between the two approaches: A\ = A l { 



hit 



Note that in the perturbative approach the sum over iV is convergent only for 
\2n\e x | < 1, but we define it in general by analytic continuation. The analytically con- 
tinued sum vanishes exponentially for large x°, which yields a convergent x° integral. 

The finite radius of convergence of the sum over N is a feature of all the amplitudes we 
have examined, and makes it challenging to extract the late time behavior. In particular, 
to find the late time behavior of the unintegrated amplitude we must first find the exact 
early time behavior, so that we have an exact formula to analytically continue. This in 
turn means that we must perform an exact summation of the perturbation series. In 
the calculation above this turned out to be easy, but the requisite sums become more 
challenging for higher point amplitudes. In light of this, it would be extremely useful to 
develop a method to access the late time behavior directly without resorting to analytic 
continuation. 

Although we have so far focused on tachyon vertex operators, it is instructive to 
consider a couple of other examples. Using the perturbative approach we can compute the 
one-point function of dX dX°e lujX at the origin to be H] 



^ e -S b ^ ydx Q dx G e iuX^ = / dx 



2 1 + 2ir\e x ° 
in (2tv\)-^ 



(3.20) 



2 sinh nu 



The S(uj) terms comes from the N = term in the sum, which is special since the bulk 
fields have nothing to contract against. On the other hand, if we repeat this calculation 
in the approach where we first integrate over x° we will miss this S(u>) term. This is the 
because the latter approach is based on analytic continuation in ui and can therefore miss 

9 This computation also applies to an arbitrary closed string state, since as pointed out in 
we can always choose a gauge such that there are no timelike oscillators. 
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non-analytic pieces. The S(u) term can be shown to be responsible for energy conserva- 
tion: it gives Tqo where T^ v is the energy-momentum tensor of the decaying brane. The 
lesson is that the integrated approach can miss certain physically relevant non-analytic 
contributions, but can still be useful, provided we keep this lesson in mind. 

Further subtleties are revealed when we study the one-point functions of higher di- 
mension operators corresponding to massive string modes. A finite number of terms in the 
sum over iV will be "special", leading to integrands which diverge for late times as e px 
for positive integer p (this behavior was observed in [^7[ , |[28|| .) Such terms will also arise 
in amplitudes involving more than one bulk tachyon, since such operators will appear in 
the OPE. This behavior clearly makes the direct evaluation of x° integrals problematic. 
To define these integrals one should impose the constraints of conformal invariance, as 
explained in [29], [|30[]. In particular, after computing the tachyon one-point function, we 
can write the remaining operators as Virasoro descendants, and then use the fact that 
boundary state is conformally invariant to generate the remaining amplitudes. 



4. Bulk-boundary two-point function: Perturbing the initial brane state 

A two-point amplitude in which one of the bulk operators is moved to the boundary 
describes an open-closed scattering amplitude. The closed string could either be an in- 
state or out-state, but because the D-brane is decaying and does not exist at late times, 
the open string vertex operator can only be describing an in-state. Thus such amplitudes 
provide a systematic way of exploring the effects of perturbations to the initial conditions 
describing a decaying brane. For example, we could construct a localized lump on the 
brane at early times and ask how long in time the effect of this deformation lasts. 

In the previous sections, we have focused on closed string tachyons. However, the 
generalization to include an arbitrary closed string state is relatively straight-forward if 
we choose a gauge in which the on-shell closed string state with non-zero energy has no 



time-like oscillators |2E[]. In this gauge, the closed string vertex operator has the form 

V = e luJcX °V sp . (4.1) 



The spatial part of this operator V sp is constructed from the 25 spatial fields and is a 

2 

Virasoro primary with conformal dimension A = 1 + The computation of a bulk- 
boundary amplitude factorizes into a product of the two point function in the time-like 
boundary Liouville theory and in the free spatial CFT. The form of the closed string 
vertex operator fl4.1| ) is useful, because in the time-like boundary Liouville CFT the vertex 
operator has the same form as the tachyon vertex operator. In the spatial direction, the 
operator is non-trivial. We can write V sp as 

Vsp = V(dX, d 2 X ■ ■ • , dX, B 2 X ■ ■ -)e l P-x , (4.2) 

for a polynomial V. In computing the bulk-boundary two point function with such an 
arbitrary closed string state, the contribution from the spatial part of the CFT is given by 

(: V(dX, d 2 X ■ ■ ■ , dX, B 2 X ■ ■ .)e^ : : e^ (t) :) , (4.3) 
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where t parameterizes the position of the open string vertex operator on the boundary. The 
integration over the zero modes in the Neumann directions will yield a factor of 5(p\\ + k\\ ). 
This 5-function will be multiplied by terms from the contractions of the bulk and boundary 
fields. For no open string vertex operator (or equivalently for k = 0) the result of these 
contractions is just a phase ||. In particular, note then that the absolute square of the 

2 

result is independent of u c , even though u> c enters into the definition of V sp via A = 1 + 
a result which is most easily seen from the boundary state. In the case with the open string 
vertex operator there will be additional contributions from the various ways of contracting 
fields from the bulk and boundary operators. This result will in general depend on the 
open string momentum k, and therefore on uj . However, the squared result will again be 
independent of u c just as before: u c dependence can potentially enter only in the Dirichlet 
part of the computatation, but since the open string vertex operator has no momentum 
components in the Dirichlet directions, this part of the computation is identical to that 
without the open string vertex operator || . So as far as the spatial CFT is concerned, the 
effect of the open string vertex operator is to contribute a ujq and V dependent factor, but 
not to alter the ui c dependence. Since our primary interest is in the u> c dependence, we will 
omit the extra uj$ and V dependence (though these could be straightforwardly computed), 
and so the formulas we will write below are strictly valid only for V = 1, i.e. for the closed 
string tachyon. 

We now return to the computation of the bulk-boundary two point function in the 
time-like boundary Liouville theory. First consider the perturbative approach. To obtain a 
correlator with one bulk and one boundary operator from ( |3.7| ), we set I — 2, use conformal 
invariance to place one operator at z\ = and the other at z 2 = 1,0 and discard the 

divergent factor (1 — Z2Z2) 2 ^ 2 since it is removed by boundary normal ordering. This 
gives 

A 2 (u; c: u; ) = J dx° e *K+^)*V(0, u c , 1, uj q \ fi), (4.4) 

with 

00 

F(0, u c , 1, u - fj)=^2 e- N »I N {u ) , (4.5) 



N=0 



and 



I N (u )= [ dU |det(l -U)\ 2luJ ° . (4.6) 

JU(N) 

Here \i is the "chemical potential" defined in (|3.9| ) and u> Cj0 are the energies of the closed 
and open string vertex operators. 



Now consider the integrated approach. Following fl3.12|) , the integrated amplitude is 



Af(o; c , uo) = -m-^rr- f , V B ^ w ») ■ ( 4 - 7 ) 

smh(7r(o; c + u ) 



10 Actually, this step is not entirely innocent. Separating out the zero mode x° breaks conformal 
invariance, and so conformal invariance can only be restored after performing the x° integration. 
Since x° is noncompact, there can be subtleties associated with boundary terms, as discussed in 



[31]. 
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Here B 2 is defined by "analytic continuation" from imaginary momenta. Setting w CiC , = 
in 0)O gives 

B 2 (in ci in ) = I N (in ) , (4.8) 

with iV = n c + n Q . Below we will evaluate A 1 ^ using both the integrated expression 
( |4.7|) and the perturbative expansion (|4.4| ) and show that they agree. This bulk-boundary 
amplitude captures the main features of more general amplitudes, but in a cleaner context. 
It allows one to vary the initial state of the brane by creating an open string perturbation. 
We can therefore explore a key physical question: are the general features of brane decay 
sensitive to the precise initial state introduced by Sen. For these reasons, we turn to a 
detailed study of the bulk-boundary two point function. 



4-1- Use of the Selberg integral 

Our basic technical goal is to evaluate integrals of the form Q4.6| ). This can be done 
with the help of a famous integral due to Selberg. (For a pedagogical review, see H32j .) 



Following Sec. 2.1 of f33 



<H7|det(l -U)\ 



-2a 



1 f dti dt 



N 



N\ J 2tt 2tt 



■>N —2aN r°° 



n 

i<j 



it ■ it 1 2 

fr 1 1 £ 3 



-it; I —2a 



N 



N\(2n) N 

The Selberg integral of interest to us is |52 



dxi ■ ■ ■ dxN J~J \xi — Xj\ 2 J^[(l + 



1=1 



(4.9) 



dxi ■ ■ - dxN ] ( 



TTfi + x 2 )-^« - {27l)Nm TT mrU ~ 2a) (4 



i<j 



1 = 1 



10) 



We thus find that the integral relevant for the open-closed amplitude is 



I N {ict) 



dU\det(l-U)\ 



-2a 



N 

n 



T(j)T(J - 2a) 
(r(j - a)) 2 



(4.11) 



Selberg's integral converges for real a when a < |, which is also the condition for the 
Gamma functions on the right side to have positive arguments. We will define the result 
for arbitrary a by analytic continuation from the convergent region. In order to define 
the integrated amplitudes (^j]) f° r g enera l rea l momenta we will have to also analytically 
continue N. The latter seems particularly problematic since the left hand side of ( 4.11 ) is 



an integral over U(N) matrices and the right side contains a discrete product involving N. 

Happily, progress can be made using the integral representation of the Gamma func- 
tion, 



iogr(z) 



00 dt 
T 



(z-l)e" 



— e 



-zt 



1-e- 



(4.12) 
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which is valid for z > 0. (This domain of validity is the same as that required for the 
convergence of the Selberg integral (|4.10|) , namely a < 1/2 for real a.) Applying this 



identity gives the result: 

r°° rjf p-Nt _ i 

log [I N (ia)) = log [B 2 (i(N - a), ia)} =J y (l - ^ ^——^ . (4.13) 

Both a and N can be analytically continued in this expression. We will use this below 
to compute and compare the perturbative (|4.5| ) and integrated (|4.7|) 2-point scattering 



amplitudes. Although the integral expression is itself only convergent in some regions of 
the complex momentum plane, we will see that it can be computed exactly in terms of 
special functions in these regions which can then be continued to general values of the 
momenta. 

4-2. Computing the Bulk-boundary amplitude 

• Integrated approach: In the approach in which the zero mode is integrated at the 
outset we wish to evaluate ( |4.7|) . We can do this by analytically continuing ( ^.13|) as 
a — > —iu and N — > —i(u c + u ). This gives: 

smh(n(uj c + uj )) 

(1 — e ~ iuj ot^2 ^ ' > 

^'^"^(l-cosht) • 

We will evaluate the integral in the exponent in closed form in terms of special functions 
below. But first let us compare Q4.14j ) to the perturbative result 



• Perturbative calculation: In the perturbative approach we must calculate the sum 
in Q4.5| ). In terms of the function H(t,u ) defined above, this becomes 



oo 

F(0, u e , 1, u 0] »)=J2 e-^e/o" * *<*.«o)(e-**-i) 

N=0 

00 p-NfJ. roc I 

j2 -n- / n dt - h ^ u °) ( e ~ Nta - v 



(4.15) 



l,N=0 " u a=l 

We can now expand out the product of (e~ taN — 1) factors in a series and explicitly carry 
out the sum on N. This gives 

1=0 m=o 11 WJo a=1 l + 2n\e x ~^ b =i tb 

(4.16) 

where we have substituted back for the chemical potential \x in terms of the coupling A and 
the zero mode x . The full amplitude (|4.4j ) involves a Fourier transform of x° as in ( f4.4|) . 
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We recognize the required transform from the expression (|3.19| ) used for the computation 
of the closed string one-point function and find 



smh(7r(w c + w )) f-^^ Q U \mj 

^° 0=1 



(4.17) 



Carrying out the sums we arrive at the final result 

A 2 {uj c ,uj ) = -in . ■ --eJo ; . (4.18) 

smh(7r(a; c + u )) 

This precisely matches the integrated amplitude (|4.14j ) For future purposes, we define 

p 00 

G{u c ,u )= dtH&UoXe^+^-l) (4.19) 
^0 

as the exponent in the bulk-boundary amplitude. 

The fact that the integrated and the perturbative approaches match exactly greatly 
increases our confidence in the methods developed here. The integrated approach is based 
on defining the expressions in ( pi.l2| ) by first expanding them for integer momenta, evaluat- 
ing the integrals and then interpreting the result as definition of the amplitude for general 
momenta. Such an "analytic continuation" is inherently problematic - without further 
specification of asymptotic conditions there exist many analytic functions which coincide 
with the one at hand at integer momenta. By contrast, the perturbative approach does not 
suffer from this failing and thus the agreement between ( [4.18| ) and ( [4.14j ) is very encourag- 
ing. Nevertheless, both approaches used here involve analytic continuation in the momenta 
and run the risk of missing non-analytic pieces. In particular, the perturbative approach 
assumed a well-defined Fourier transform and thus could miss non-analytic contributions 
to the amplitude associated with zero momentum processes. Despite this we will be able 
to extract useful insights into the structure of D-brane decay from ( ^.18|) . 

4-3. Comparison with other approaches 

As discussed in the previous section, one might have proposed an alternative approach 
to evaluating these amplitudes, by analytically continuing 0J2 to imaginary integer values 
(—iiV2 = n G Z + ), since in that case jl~6|1 shows that 



f dU \ det(l - zU)\- 2n = (1- zz)~ n2 . (4.20) 

JU(N) 



11 Related results in the Liouville approach appear in [IS]. 



14 



This expression is valid for n and N nonnegative integers with n < N. Except for the 
trivial case iV = 0, the region where this formula is valid is disjoint from the region a < \ 
where the Selberg integral was valid prior to analytic continuation in a. To compare the 
two approaches, set a = n G Z + in Q4.11|) and take z = 1 in ( |4.20|) to get the proposed 
answer for the bulk-boundary amplitude. We find the following two possible results 

n™-?= finite, 
dU | det(l - U)\~ 2n = { f = \ TO - n)) 1 (4.21) 



(0)-" =oo . 

The first expression is obtained from the Selberg formula used in (^4.11] ) while the second 
arises from ( |4.2C| ) as derived in The Selberg result is finite for N > 2n and vanishes 
for iV < 2n. By contrast, the second formula diverges for any N. The consequence of this 
is that after continuation back to real momenta, naive use of the latter expression from 
||16|| leads to a vanishing amplitude for all momenta, while the Selberg formula will lead to 
a finite result. 

This situation can be given a useful worldsheet interpretation. When a bulk vertex 
operator is taken towards the boundary to define a bulk-boundary amplitude, it will collide 
with operators from the boundary interaction. Therefore, we need to add counterterms 
to operators which are taken to the boundary in order to dress them appropriately in 
the presence of the interaction. The second formula in ( |4.21| ) corresponds, however, to 
a bare vertex operator, and so diverges. On the other hand, the analytic continuation 
used to define the Selberg integral is a convenient way of regulating and renormalizing 
the operator, and so gives a finite result. Of the two results displayed in ( fl.21| ), the finite 
Selberg result is the physically relevant one. 

The summation over N that was carried out in the perturbative approach, and which 
led to ( f4.18| ) for the final amplitude, became possible because the Fourier transform in 
the definition of A2 Q4.4|) simplified the expressions. However, it can also be of interest 
to do the sum on N in ( |4.5| ) without integrating over the zero mode a; . As we discussed 
earlier, the contributions to the scattering amplitude from late times (large x°) could shed 
light on the nature of the "tachyon matter" state to which the brane is supposed to decay. 
Doing the perturbative sum in ( |4.5| ) exactly for fixed small x° and then continuing to large 
x° turns out to be very difficult for arbitrary momenta. However, it is possible to make 
progress for special imaginary integer momenta. The relevant techniques and results are 
presented in Appendix A. 

4-4- Initial state perturbations and the decay of the brane 

Our result ( |4.14| ) (or ( f4.18| )) for closed string emission from a decaying brane perturbed 
by an additional open string tachyon allows us to explore how the emission of closed strings 
is modified when the initial conditions for the brane are changed. For example, we might 
make a tachyon lump on the brane by superposing open string vertex operators and ask 
how this affects the decay. 

Important questions to answer are whether all the energy in the brane decays into 
closed strings, and what is the distribution of these strings. Sen has argued that brane 
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decay leads to a pressureless state of "tachyon matter" 0] . A computation of closed string 
emission from decaying branes by || (summarized in eq. ( |3.19| ) in Sec. 3) showed that at 
tree level the emission of closed strings is exponentially suppressed in their energy. For 
low dimension branes this nevertheless implies a divergent total emission, because of the 
Hagedorn growth in the density of states. This divergence was interpreted as indicating 
that all the energy of the decaying brane would be converted into closed strings, and that 
very heavy strings would dominate the decay products. Such heavy strings would have 
the stress tensor of pressureless dust, thus suggesting that they constitute the mysterious 
tachyon matter. It was also argued that near the endpoint of the decay, including back- 
reaction would ensure that energy conservation was satisfied so that the emission of closed 
strings shuts off after all the energy in the original brane has been converted. This was 
confirmed in the c = 1 matrix model, where it was also found that it is crucial to take into 



account the quantum mechanical nature of the decaying brane | 12]| . By contrast higher 
dimensional branes had a finite total emission of energy. The conventional wisdom states 
that in this case, small perturbations can also lead to decay into higher co-dimension branes 
which would account for the "missing energy" in the decay into closed strings. 

Our results can enable a systematic exploration of how changing the brane initial 
conditions affects the decay products. To initiate this study we will compute a closed form 
expression for the closed string emission amplitude ( 4.14Q in terms of special functions and 



then extract the asymptotics of the decay for large closed string energies. 

Before plunging into the calculation, let us note that there are two possibilities for what 
one might mean by "changing the brane initial conditions". One interpretation involves 
describing the brane by a boundary state; then to change the boundary state we can act 
with an exponentiated open string vertex operator (at least to first order in the operator; 
beyond that there will be corrections). This corresponds to changing the classical open 
string tachyon profile of the decaying brane. To first order in the perturbation, the total 
emission rate into closed strings will then be the sum of the original rate plus a correction 
due to the perturbation. Here we will focus on an alternative interpretation, in which 
we add a single quantum open string excitation on top of the classical background. The 
open-closed amplitude then corresponds to the amplitude for the incoming open string to 
disappear and for a closed string to be created. 

To compute the bulk-boundary amplitude ( f4. 14| ) in closed form we must evaluate an 
integral of the form 

G M =l - il -e-^ Ir —^ ; /J=-i (ae+ «.). (4.22) 

Notice first that this integral vanishes if uj = 0, i.e., in the absence of any open string 
perturbation. In that case, 

A 2 {o Jc ^) = -m \ , (4.23) 

smh(iT(jj c ) 

reproducing the known result in ||. So the full result has the form of a modulation of the 
emission amplitude without an open string perturbation by a factor of e G . It is also easy 



16 



to argue that the real part of G will be negative so that the modulating factor suppresses 
the amplitude. 

The challenge in evaluating the integral ( f4.22| ) occurs because as t — > 0, cancellations 
between the terms in the numerator are needed to cancel the overall power of t 3 in the 
denominator in this limit, leading to a finite integral. With this in mind, to evaluate G for 
general momenta, observe that G also vanishes if (3 = i(u c + iv ) = 0. Because of this, we 
can write 



d(3 



G(u c ,u )= d/3— =- dp -{l-e-^Y- — . (4.24) 



o 2 v ' (1-coshi) 



After expanding (1 — e luJot ) 2 , each term in dG/df3 can be evaluated with the help of the 
integral formula 



Jo \ > r( 7 + , + i) 2 



Re, > -- , (4.25) 

£4 

which is listed in |34]] (vol. 1, p. 163). This expression diverges as v —* —1, but the linear 
combination 



■ i ~ i ' , Lr/0l , .. L | J r (7 - v) _ 2 rh-v + iu ) r( 7 - v + 2iuj ) 

T(7 + i/ + l) r( 7 + i/ + 1 + zw ) r( 7 + v + 1 + 2iu a ) 



(4.26) 

needed here is well-defined as , — > —1. 

The limit v — > — 1 gives dG/d(3 in terms of derivatives of the Gamma function. 
The integral with respect to j3 that remains can be evaluated by again using the integral 
representation ( f4.13| ). The result can be written in terms of the Hurwitz £ function 

1 fOG j.s-1 -zt °° 1 

r(a) Jo 1 - e * ^ (n + z) s 

and £ m ' n (s, z) = g s mg z n ■ In order to summarize the results, we introduce the notation 

H[F(x), a, b)] = F{b + a) - 2F(b) + F(b - a) . (4.28) 
Then, after using various £ function identities, we find the result: 

G{u c ,u ) = W[C (1 ' 0) (-1, x), -iu , iu ] - H[C (1 ' 0) (-1, x), -iu , -iu c ] 

+ H[x(((0, x) — lnT(x)), — iw , iu ] — TC[x(£(0, x) — lnr(x)), — iu , — iu c ] . 

(4.29) 

Putting this expression back into ( f4.18| ) gives the general bulk-boundary amplitude in 
closed form in terms of zeta functions. This result can now used to extract many aspects 
of the physics of decaying branes. Here we will explore one question - how does including 
an open string perturbation change the asymptotics of brane decay into closed strings? 
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To do this we must find the asymptotics of G(u c ,u! ) for u c ^> u a . We will use the 
expressions 

C(o,x) = i -x , 

(10) B 2 (a)\na B 2 (a) + a B 2 ^ a~ k + 2 B k ( 4 ' 3 °) 

C 1 ,GJ 2 4 + 4 ^ (fc-2)(fc-l)fc ' 

Here B 2 (a) = a 2 — a + 1/6 and I?2 = 1/6 are the second Bernoulli polynomial and number 
respectively. The second formula was derived for 3?(a) — > oo fl3~5]l , but it will apply by 
analytic continuation for large |a| in general. Applying these formulae we find that for 
uj c 3> w , the dominant term is G(u> c , uj ) ~ — u 2 lnw c . If u is itself large the formula is 

G(^ )^-^ 2 ln(^) . (4.31) 

For small u a and u c 3> u; the denominator in the log is modified, but the —U3 2 \nu3 c 
behaviour always holds. The next to leading terms are independent of lo c , but dependent 
on u> and include a phase. Beyond that terms are suppressed by powers of u c . We have 
checked these asymptotics numerically. 

To summarize, the amplitude for emitting high energy closed strings, when the brane 
initial state has been perturbed by a reasonably energetic open string which is subsequently 
absorbed, is 

A 2 (u c ,u ) ps -in ) / ■ rr — ■■• (4.32) 



sinh(7r(c<; c + u a )) \uj. 

where the ellipses represent terms that are either constant or fall off as uj c increases. 

For the unperturbed initial state of the D-brane (the spatially homogeneous decay), 
it was found in || that the total energy of closed strings emitted was divergent for a Dp 
brane with p < 2 (but see also [^6] ). We now examine how this statement is modified 
when the initial state is perturbed by addition of the boundary tachyon vertex operator. 
To compute the expectation value of total emitted energy we add up the squares of the 
individual emission amplitudes, 



^~E|l4 S Wo)| 2 ■ (4.33) 

For large u c 



\A 2 (u c ,uj )\ ~ e - 2 ™°u-» 2 ° . (4.34) 

—uj 2 

This differs from the unperturbed case in the extra factor of uj c °. A similar calculation 
to the one done in H yields 



E 1 

Vp ~ W) 



(4.35) 
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Notice that the extra factor of uj c ° is a suppression if uj q is real but an enhancement for 
imaginary uj . Without the initial state perturbation, the average energy per unit volume 
diverges for the DO, Dl, and D2 branes. For the DO brane, the on shell condition fixes 
the initial state perturbation with the open-string tachyon vertex operator to be e x , i.e. 
uj = i and the perturbation enhances the closed string emission leading to a divergent 
result. For the Dl brane, the result depends on the details of the perturbation. If the spatial 
momentum is large enough, such that ujq > 1, the closed string emission will be finite. For 
higher dimensional branes, the total emission will be finite for any perturbation with real 
ujq. If the spatial momentum on any brane is sufficiently low, the on-shell condition will 
cause uj to be imaginary, and there will be an enhancement of closed string emission. For 
p < 3 the energy contained in closed string emission will diverge in this case. 



4-5. Comments 

For the DO brane, the divergent closed string emission is expected to be an artifact 
of perturbation theory, which will be regulated by back reaction and quantum effects. It 
strongly suggests that the initial D-brane decays completely into closed strings and that 
tachyon matter, the end product of the decay in classical open string theory, is made up 
of very heavy closed strings |§ . A similar interpretation is given to the divergent emission 
of closed strings from the unperturbed Dl and D2 branes. For higher dimensional branes 
(p > 3) the conventional wisdom has been that the decay process is more elaborate since it 



proceeds inhomogeneously (0, |37[| , flog]). Notice, however, that in a quantum mechanical 
treatment of a decaying brane of any dimension, there will always be a component of the 
brane wavefunction along operators with imaginary ui Q . This is because such an open 
string operator will be proportional to e^°^ x and so be very small at early times, in effect 
corresponding to an infinitesimal perturbation of the initial state. We have shown that 
such perturbations enhance closed string emission for all Dp-branes, even when the decay 
is homogeneous, and would lead to divergent emission (without treating backreaction and 
wavefunction effects) for p < 3. In fact, we will get a divergent emission for any p if we 
are willing to consider non-normalizable perturbations with imaginary spatial momenta. 
This enhancement of closed string emission is appealing. It suggests that a full treatment 
of decaying branes of any dimension will show all the energy of the brane entering into a 
state of heavy closed strings. 

In our study of the open-closed amplitude we have discussed the open string state as an 
in-state, or a perturbation of the initial condition, while the closed string is emitted as an 
out-state. By complex conjugating we could have naively written down an amplitude that 
would apparently describe the emission of an open string. Likewise, at the perturbative 
level, we could have computed amplitudes containing additional strings in the final state 
including, apparently, open strings. This seems puzzling given that the brane has decayed 
at late times and thus open strings should not exist. It is possible that such amplitudes are 
simply inconsistent at a non-perturbative level. However, it seems more reasonable that we 
should interpret such amplitudes with "open strings" in the final state as follows. In any 
time dependent background we construct vertex operators at early and late times (in- and 
out- operators) to describe simple perturbations around the respective vacua. However, we 
are always free to use the complete basis of in-operators to describe outgoing states also 
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- these operators will simply describe very complicated out-states. In our context, both 
open and closed strings can exist at early times, but only closed strings exist at late times. 
An out-open-string operator should thus represent some very complicated correlated state 
of closed strings that emerges from the open string via the decay of the brane. This 
interpretation leads to the important question of whether, in order to compute the full 
decay of a brane, we must calculate amplitudes with all possible out vertex operators, or 
whether computing the emission of closed strings as in ||, or as done above, suffices. In 
other words, do we include all decay channels that are perturbatively apparent including 
the ones with final-state "open strings" , or is this "double counting" ? The answer to this 
question goes to the heart of the open-closed duality for decaying branes that has been 
conjectured by Sen (JT^j, 0). Regardless, the amplitudes computed above are essential 
components of the full answer and are the complete result for the "exclusive" amplitude 
describing absorption of the initial open string and emission of the final closed string. 



5. The bulk-bulk amplitude 

Now we turn to a computation of the bulk-bulk amplitude which describes the scatter- 
ing of a closed string from the decaying brane, or the emission/absorption of a correlated 
pair of closed strings during the decay. To obtain the relevant correlator from (|3.7|) and 
( P.8|) we use the conformal symmetry (but recall footnote 10) to locate one vertex operator 
at z = and the other at z = r, < r < 1 on the unit disk. This gives 

A 2 (u u u 2 ,r) = |r|~™ (1 - |r| 2 )-^ 2 / 2 J dx° e lx ° ^ 1+UJ2) F(0, u lt r, u 2 - n) , (5.1) 

with \i = —x° — ln(— 2tt\) and 

oo 

F(0,u 1 ,r,u 2 ;iJ,) = e~ NfJ, J N (u 2 ,r) , 

N=0 

J N (u 2 ,r)= I dU \ det(l - rU)] 2 ^ 2 . 

JU(N) 

To compute the full amplitude we should integrate over r as J^drr. 

5.1. The Toeplitz determinant 

The key steps in computing the amplitude are to evaluate the terms Jn in the infinite 



(5.2) 



series (|5.2| ) and carry out the sum. Each Jn is an expectation value of a periodic function 



with respect to the circular unitary ensemble (Appendix B): 



N - ° =1 a< " (5.3) 



Eat 



20 



By Heine's identity [39|, this expectation value is equal to the Toeplitz determinant of the 
Fourier coefficients of /, 



n 

E*{n /(*<)} = D *lf] = dfit(/i-fc)i<j,fc<JV • (5-4) 
i=i 

By this notation we mean the determinant of a matrix in which the entry in the j th row 
and k th column is the (j — k) th Fourier coefficient of /. For reference, the proof of this 
identity is included in the Appendix B. In our case the Fourier coefficients are (denoting 
a = —iu)2) 

or _ re _ lt 2ae _ lkt = + ^ + r2 

Z7r \k\B(a, \k\j 

Here B is a beta function and F is a hypergeometric function. This integral representation 
of the hypergeometric function is valid for a / 0, —1, —2 • ■ ■ and \r\ < 1 [ fTOfl . As we will 
see below, the boundary limit r — > 1 of the right hand side will be well-defined anyway. 
The restriction on a can also be written as uj 2 ^ 0, — i, —2i, —3i 

• Special imaginary momenta: When c^2 = ^, 2z, 3z • • • (i.e., positive imaginary integer 
momenta), or equivalently a = 1, 2, 3 • ■ • with a < N we know that the integral ( p.2| ) could 
have been done by the group theoretic methods of [If]] with the result that 

J N (uj 2 ,r) = (1 -r 2 ) u i ; w 2 = in, n G ; n<N. (5.6) 

Since these momenta are within the range of validity of the integral representation (]5.5|) , 
the determinant Ujvt/] i n ( |5.4| ) must reproduce ( |5.6| ) at positive imaginary integer u>2- It 
is easy to verify that this is indeed the case, at least for small n and N but we have not 
constructed a general proof. 

• Boundary limit gives Selberg: As r 2 — > 1, the integrals and coefficients fl5.4j ) reduce 
to ( f4.11|) of the previous bulk-boundary calculation. In the limit, the determinants ( |5.4| ) 
reduce to 

/ r(l-2a) \* ^ r(a+\i-j\) 1 
Hence this determinant must be equal to fl4.11| ): 

j N {u 2 , i) = D N [f\ r=1 = f[ r y r(j ~ y . (5.8) 

7 = i ( r u - «)) 

We have checked the identity explicitly with Maple for = 1, . . . , 8, but have not con- 
structed a general proof. 

These two special limits appear initially to be in contradiction with each other since, 
for u>2 = in, (|5.6| ) diverges as r — > 1, while ( |5.8j ) is finite. (Also see the discussion around 
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( [4.21 ).) However, what we really seeing here is that the order of limits can matter when 



one approaches the boundary at r — > 1 at special momenta. Our Toeplitz determinant 
( p.4| ) provides a definition of the two point amplitude for general momenta and gives a 
well-defined analytic continuation to the bulk-boundary amplitudes (with r = 1) that 
we derived in the previous section. In addition the determinant reproduces the special 
form (|5.6|) which was previously obtained at imaginary integer moment |TH]. Therefore, 



it appears that naive analytic continuation of ( |5.6| ) to general cu2, which is the technique 
used in the past, is invalid, as one might have expected. 

Below we will see that the naive continuation of ( |5.6| ) from imaginary integers to 
general uj2 amounts to a large N approximation. 

5.2. Large N approximation 

In general, dealing with the Toeplitz determinants is very complicated. We will con- 
sider two methods. For large N, approximate results can be found by using Szego's limit 
theorem. This estimate can be improved by using an identity due to Borodin and Ok- 
ounkov HU (see also P2||). The latter authors considered Toeplitz determinants D n [f] 
of periodic functions and related them to Fredholm determinants of certain kernels. In 
particular, they considered functions of type 

/ = (i+£ l77 )*(i+£ 27r y , (5.9) 

where £i, £2) z -, z ' are complex parameters and r\ is a coordinate on a complex plane. If we 
set 77 = — e -1 *, £1 = £2 = r, and z = z' = —a, we recover our particular case. Here a can 



be an arbitrary complex number. Following [41], there is an identity 



D N [f]=Z det(l-K N ) . (5.10) 

The prefactor Z involves the Fourier coefficients of In/: 

00 

Z = exp(£fc(h/) fc (£/)_ fc ) . (5.11) 
k=i 

This is the large iV approximation for the determinant given by Szego's limit theorem, 
lini/v^oo Djtf[f] = Z, whereas the factor det(l — Kn) contains all the corrections to the 
approximation. To evaluate Z we use the relation [f40| (assuming < r 2 < 1) 

ln(l - 2r cos(t) + r 2 )e- ikt dt = -a—, (5.12) 

-7T 1^1 



which gives 

D N [f] w Z = exp{a 2 ^ k r 2k k~ 2 } = (1 - r 2 )"° 2 . (5.13) 

k=i 

This large iV approximation reproduces the naive analytic continuation n 2 — > a 2 of ( |3.16| ) 



which results from calculations at imaginary integer momenta [16]. In fact, as we showed 
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in (|5.6|), when a = n, n G Z + , n < N, the result ( |5.13| ) is exact, or equivalently, the 
Borodin-Okounkov determinant det(l — K^) equals 1. There appears to be a remarkable 
localization phenomenon in the amplitude ( |5.2| ) such that for special momenta the result 
localizes to the large N result Z (|5.11| ). 

Notice that the right hand side of (|5.13|) is independent of the dimension N of the 



matrix. It is dangerous to use this approximate large iV result in the bulk-bulk amplitude 
(pT2|) because of the infinite sum, but if we do so regardless, we get (restoring a = —iujq) 

OO /.. 2\lu 2 

F(0,c*,r,a»;/i) = £ e"^(l - r 2 )^ = { ' • (5-14) 

N=0 



Putting this back in the complete amplitude (5.1) and integrating over x° gives 



A 2 (u 1 ,u 2 ,r) = -m4^ 7 ■ -\r\-^(l -\r\ 2 )^' 2 . (5.15) 

sinh(7r(o;i + uj 2 )) 

This gives the same result as in (|3.18| ), computed by doing the zero mode integral first 
and using naive continuation from large imaginary integer momenta. The full amplitude is 
computed by including the spatial part of the amplitude and integrating over the position 
of the vertex operator. To do this recall that the spatial part of the conformal field 
is unchanged by the boundary perturbation, so we can include the contributions for a 
standard D-brane in ( |2.11| ), and write 

A 2 = / drrA 2 (u 1 ,u 2 ,r) = -m \ - >, -/ r IrfM - H 2 ) s/2 ~ 2+ " 2 , (5.16) 

Jo ~ smh(7r(o;i + u 2 )) J 



where we used the Mandelstam variables ( |2.2p and the closed string tachyon on-shell con- 
dition (( |2.3| ) with N = 0). Finally, integrating over r gives the result 

. (27rA)-^+^) r((t/2) + i)r(£ 2 -i) 

A 2 = -iTT . \ . '-, ^ , (5.17) 

smh(7r(^i +u 2 )) T((t/2) + fc 2 ) 

where k\\ is the spacelike part of the momentum in the directions parallel to the brane. The 
factor r(£/2 + 1) in this expression contains the closed string poles ( |2.3|) . The remaining 
poles coming from the T(/c 2 | — 1) are more exotic, since they occur at special spacelike 
momenta. If momentum here was not simply the spacelike piece, we would be recovering 
open string poles. 

Above we pointed out that the boundary limit r — > 1 of the full amplitude (|5.8| ) gives 
the result of the Selberg integral that we used to compute the bulk-boundary amplitudes in 
Sec. 4. By contrast, the large N approximation described here, since it coincides with the 
naive continuation from imaginary integer momenta, gives a singular boundary limit. (See 
Q4.21 ) and the discussion around it.) In Sec. 4.3 we argued that the singular behavior of the 



naive analytic continuation could be related to a need to regulate and renormalize bulk 
operators as they approach the boundary and collide with the boundary perturbation. 
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Here we have seen that a large iV approximation of the bulk result leads to a singular 
boundary limit, while the full expression leads to the finite quantity that was used in 
Sec. 4 to study the effect of initial state perturbations on the decay of the brane. It may 
be useful to think of the large iV approximation ( |5.13|) as a saddlepoint contribution, with 
the Borodin-Okounkov determinant in (|5.10|) containing corrections due to fluctuations in 



the interactions between the vertex operators and the boundary perturbations. Perhaps 
it is possible to understand the determinant in terms of a renormalization or dressing of 
string vertex operators by the boundary interaction. 

Our result demonstrates that the usual technique for computing decaying brane ampli- 
tudes by analytically continuing from imaginary integer momenta is unreliable. The naive 
continuation only isolates a large N contribution to the full answer. This happens because 
of the remarkable localization of the full integral to the large N contribution at imaginary 
integer momenta. Our amplitudes are directly defined at real momenta and differ in im- 
portant ways from the naive continuation. We have already discussed the consequences of 
the extra contributions for open-closed amplitudes. We now turn to a brief examination 
of the extra factor, the Borodin-Okounkov determinant, in a general amplitude. 

5.3. The Borodin-Okounkov determinant 



The determinant in the full amplitude Disr[f] = Z det(l — Kn) can be expressed ( [fiT 



and ||42|| ) as a Fredholm determinant of an operator K^, acting on the Z 2 normed infinite 
vector space Z 2 (iV, N + 1, . . .), defined by 

00 00 

det(i -k) = i + j2 (-ir E de W- ■ ( 5 - 18 ) 

m=l N<h<l 2 <---<l m 

This equation is essentially a discrete version of the Fredholm determinant of a kernel 
familiar from the theory of integral equations. For the particular case of ( |5.9|) (with the 
parameters replaced by ours), the kernel K takes the form (correcting a minor error in 

H) 



K(i,j) = (-^W-a) J+V +j+2(1 _ 7 , 2) - 2q _ 1 
i.j . 

1 . . . r 2 ,F(-a- 1,1 + a,j + 2; 7 £_ 1 



x 



(F(a, -a, i + 1, ^— r)- '— — p (5.19) 

r z — 1 7 + 1 



i — j r 2 — 1 j + 1 

F(-l - a, 1 + a, i + 2 



2 



i + 1 



F(a,-a,j + 1;^— -)) 



The notation (0)5 indicates Pochhammer symbols. From the point of view of the scattering 
amplitude, an important property of the kernel K is that it contains the information that 
is missed in the naive analytic continuation of the previous results. 

As a consistency check, we can examine the behavior of (|5.18|) as r — > 0. From 



we can see that the Toeplitz determinants reduce to 1 as r —>■ 0. In addition, from ( |5.13| ) 
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Z — ► 1 as r — > 1, so we expect that det(l — K) -> 1 as well. In the limit the kernel 
factorizes: 

KikJj) ^ r 2 ■ (-a) h+1 r h (-a) h+1 r l i . (5.20) 

As r — > this of course implies that det(l — K) — * 1. In fact, a stronger statement is true 
- det K vanishes to leading order in r: 

m 

&et[K{l^ 3 )}™ =1 = r 2m \[{-a) h ^ 

4=1 (5.21) 
x (-in-a) K{1 /^ ■ ■ ■ {-a) Um /^) = o . 

Because of the this the det(l —K) factor has no effect on the pole structure of the amplitude 
at r = 0, in agreement with what we expect from the OPE of the tachyon vertex operators. 

We already know that in the r — > 1 limit the Toeplitz determinant reduces to the 
form (|5.8|) and the bulk-bulk amplitude reduces to the bulk-boundary amplitude. Since 
the expression Z ( |5.13| ) is singular as r — > 1, there must be a compensating singularity in 
det(l — K). This is more challenging to isolate. It would be very interesting to understand 
the properties of the determinant better since many aspects of the physics are clearly 
stored in it. 



5.4- What we learn 

We have given a prescription for defining the two point scattering amplitude for decay- 
ing branes at general momenta, and demonstrated that the usual tactic of naively contin- 
uing from imaginary integer momenta is a large N approximation that misses important 
contributions. The naive continuation fails because, remarkably, the general amplitude 
localizes to the large N contribution at these special momenta. Since we are using ana- 
lytic continuation in our methods it is also possible that non-analytic contributions such 
as zero-momentum delta functions are missed. These are also essential to understand, 
because they make contributions to important quantities such as the stress tensor. In any 
case, it is clear that understanding the properties of Toeplitz determinants emerging from 
U(N) random matrix theory will lead to progress in the study of decaying branes and vice 
versa. 

A useful way of making progress with the bulk-boundary amplitude in Sec. 4 was to 
do the integral over a; first, thus avoiding the need to carry out a perturbative sum over 



N as in (|5.2| ). Instead, the result then depends on a single U(N) integral with iV being 
related to the total energy in the process. Such an approach might be useful here since for 
large energies the simple expression ( p.!3| ) might suffice. 



The general amplitude prior to integrating over the time x° is: 

oo 

F(0,u} 1 ,r,u 2 ;fj,)=J2 e ~ NfMJ N(^2,r) ; \i = -x° - ln(-27rA) . (5.22) 



N=0 



This has the form of a computation in a grand canonical ensemble of U(N) matrix models. 
At very early times (x° — > — oo) the expression is dominated with smallest matrices. As 
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a; increases, larger and larger matrices begin to play an important role. As a; increases 
beyond — ln(27r A) infinitely big matrices dominate the result and the formal power series 
fails to converge. To find the contribution to the amplitude from such late times, we must 
analytically continue the fully summed expression from convergent region x° < — ln(27rA). 
In view of the many relationships between large U(N) matrices and closed strings, as well 
as recent developments in 2d string theory [fT2f , fHf ], flnfl), it is tempting to conclude 

that as the brane decays and spacetime becomes more closed-string-like in character, large 
U(N) matrices emerge from the dynamics of the disappearing open strings to describe the 
physical degrees of freedom of spacetime. 



6. Discussion 

Several interesting questions in the physics of decaying branes, including the nature of 
the final "tachyon matter" state and the possibility of defining a decoupling limit leading 
to timelike holography, require a deeper understanding of how strings scatter from such 
unstable objects. This has been a difficult problem to solve because the decay involves a 
non-trivial boundary interaction on the open string worldsheet. The standard technique 
to calculate amplitudes in this system is to first perform the necessary integrals at a spe- 
cial discrete set of momenta and then optimisitically "analytically continue" the resulting 
formal expressions. Of course this procedure is not strictly well-defined - an analytic 
function cannot be defined from its values at a discrete set of arguments without further 
information, such as either the behavior at infinity or consistency conditions derived from 
symmetry arguments. However, in the absence of further information this is the best we 
can do. In this paper we have made progress in arriving at a better-defined prescription 
for computing amplitudes. 

First we showed that the timelike part of general amplitudes could be written in terms 
of correlators in U (N) matrix models. (While we only studied tachyon scattering, vertex 
operators with more oscillators could be treated in a similar manner.) In a perturbative 
approach, a grand canonical ensemble of matrix models appeared, while integrating over 
the zero mode led to a single matrix expectation value. A key point is that the resulting 
open-closed 2-point amplitude can be exactly computed in an open subset of the complex 
momentum plane using Selberg's integral. Analytic continuation from this region to real 
momenta is reliable. Indeed, we were able to compute the open-closed amplitude in closed 
form in terms of zeta functions, and we derived the asymptotics of closed string emission 
when the initial condition for the decay is perturbed by an open string operator. In par- 
ticular we found that for extended branes the emission of heavy closed strings could be 
enhanced or suppressed by the presence of the additional vertex operator, with interesting 
consequences for the structure of the "tachyon matter" state to which the brane decays. 
We likewise showed that the closed-closed scattering amplitude can be expressed exactly 
in terms of certain Toeplitz determinants of hypergeometric functions. The results ob- 
tained by naive analytic continuation from integer momenta amount to the leading large 
N approximation (due to Szego) of these determinants. In time dependent quantum field 
theories, different choices of analytic continuation in amplitudes are often related to differ- 
ent choices of vacuum state. We did not address the issue of the vacuum in detail here, but 
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it would interesting to relate choices of analytic continuation to choices in the definition 
of the boundary state for the decaying brane. (See [§, |43|1 , |19|| , [pQ| , for some relevant 
discussion.) 

In this paper we have focussed on the bulk-bulk and bulk-boundary two point func- 
tions, and it is an important problem to study more general correlation functions. For 
example, it would be interesting to study the boundary-boundary correlation function in 
depth using the ideas developed here. One would like to make contact with the Liouville 
based approach in, e.g., (H, p|l, pfl, fl20fl, plj). 

In addition to the techniques described in the main text there are other powerful tools 
that can be brought to bear on the subject of scattering from decaying branes. In Ap- 
pendix C we show, by extending old techniques of Douglas E3] , that there is a free fermion 
formulation of scattering from decaying branes. This is related to the free fermion tech- 
niques that can be used to study decaying branes in 2d string theory [fl2| , flj~3|] , Hl4j] ). 
Appendix D points out the relation between our scattering problem and a grand canonical 
ensemble for a classical Coulomb gas. It would be very interesting to apply these methods, 
and the techniques we have developed in the main text, to a study of tachyon matter 
and to explore potential decoupling limits leading to timelike holography, perhaps in the 
context of recently discovered non-singular spacetime backgrounds representing decaying 
branes (@, |5[, @ 
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Appendix A. Special momenta and summing over N in the bulk-boundary 
amplitude 

At late (early) times x°, the chemical potential \i ( |3.9| ), is large and negative (positive). 
Thus at late times, the series expression (4J3) combined with the Selberg result Q4.11Q for 



the bulk-boundary amplitude does not converge. To get an explicit function of a; we must 
sum ( [4.5D at early times and then analytically continue. Such a procedure is rather delicate 
- terms that are small at early times may dominate after the analytic continuation to late 
times, so the sum needs to be done exactly, or at least with attention paid to the pieces 
that dominate late time behavior. Progress can be made for special momenta. 
We can write the amplitude ( (Of ) as 

f(o, a, 1? i, u 2; A = i + f; x N n r ffr + + 2 ^ 2) = 1 + E * N < N > «*) > (A- 1 ) 

N=l i = l ^ U+^JJ N=1 
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with x = e M = — 2Tr\e x . The coefficients obey 



c(N + l,u 2 ) _ T(N + 1)T(N + 1 + 2iiv 2 ) 
c(N,u 2 ) 



{T(N + 1 + iu 2 )Y 



(A.2) 



A series of this form is, by definition, a hypergeometric function if the ratio ( |A.2| ) is a 
rational function of N. In our case, this happens if iu 2 is a positive integer. In fact, when 
iu 2 = m is a non-negative integer a little algebra shows that c(N, —im) is a polynomial in 
N: 



N 



c(N, -im) = Yl 



T(j)T(j + 2m) 
y (T(j + m)) 2 



1 + 



N 



m 



m—l r 



n 

k=i 



(1 + ¥ )(1 + 2^ 



(A.3) 



Furthermore c(0, -im) = 1. 

Setting iu 2 = m the amplitude becomes — im; fi) = 1 + 5^jv = i x N c(N, — im) 
which we can be written 



d 



F(0, ui, 1, u 2 ; fi) = 1 + c ( x ^' -im) ^ 



x 



AT 



AT=1 



1 + c(x — , —im) . 

ax 1 — x 



(A.4) 



We thus arrive at 



F(0,wi,l,-zm;^) = 1+1 + 



x 



m m-l 



m 



n 

k=l 



-I k 



x 



(1 + -^)(1 + 



X 



dx 



2m — k' 



x 



1 — x 



(A.5) 



The simplest cases are 



F(0,wi,0,0;/i) 
F(0,wi,l,-z;^) 

F(0,wi,l,-2i;|i) 



1 



1 — x 
1 



(1-x) 2 ' 

1 x 3 — 5x 2 + x — 3 
~3 (1 -x) 5 



(A.6) 



Note that the late time region is given by lim 



xq— >oo 



X 



that 



lim F(0, u>i, 1, — im; fi) = 0. 

XQ — >oo 



-oo. It follows from (|A.5| ) 

(A.7) 



So we find that for these special values of momenta the bulk-boundary amplitude vanishes 
at late time. 

To summarize, ( |A.5|) combined with (|4.4|) gives the result for the bulk-boundary two- 
point function at special momenta. 
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Appendix B. Heine's Identity 

We consider integrals of the form 

i=l i=l i<j 

containing a 2-k— periodic function f(t). In the Random matrix literature this corresponds 
to an expectation value of a function Y\7=i /(^») with respect to the circular unitary en- 
semble P2|. In our case, fit) = |1 — ze~ zt \ 2,lk . A well known identity due to Heine 
|3~9[j , states that the expectation value of a periodic function fit) with respect to the 
circular unitary ensemble is equal to the Toeplitz determinant of the Fourier coefficients 
fk = J f(t) exp(-ikt)dt, 

^ = D n [f] = det(/ (i _ fc) ) 1 < ijfc < n . (B.2) 

A Toeplitz determinant is a determinant of a (Toeplitz) matrix where all entries are equal 
along diagonals. For convenience, we review the proof of Heine's identity ||39|| . 

First, note that the rii<j \ elti ~ e%tj | 2 term is equal to AA where A is the Vandermonde 
determinant A = detle 1 ^^ }k,i=i,---n- The integral can be written as: 

•' i=i 1=1 

Now expand the determinants out explicitly, A = Xln( — l) n e l ^-' i = 1 n< -^ t! , where II is a 
permutation of integers 1 ■ ■ ■ n. Hence, 

AA = £ n 1 (-l) n .+ n *e'( (n ' (,, - na(,,) ") ■ (B.4) 
ni,n 2 1=0 

Inserting this in (|B.3|), we obtain 



In=J2 ("l) ni+n2 fl { / Y~ m e l[(ni( °- n2(0)t 4 • (B.5) 
rii,n 2 1=1 ^ J ^ > 



Let vi = Ui(l) - n 2 (/). Then 



I n =J2 (-l) ni+na fl { / p-f(ti) e*»*« ) . (B.6) 
n!,n 2 i=i ^ J 71 > 

We can fix the first permutation Hi = 1, 2, 3, • • • , n and multiply the result by nl. This 
simplifies the above expression: 



n , 

I. = n!^(-l) n n /(«-n(0) 
n i=i ^ 

= n\ det(/(;_ m ))!<; jm < n . 
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(B.7) 



The determinant in the final expression is a Toeplitz determinant of the Fourier coefficients 
af/(t). 



Appendix C. Fermionization 

In this appendix, we outline a relation between the closed string scattering amplitudes 
and correlation functions in a theory of free fermions on a circle which gives an interesting 
reformulation of the scattering problem. The basic idea is due to Douglas [23] who showed 



the relation between U(N) group theory and a theory of N free fermions on a circle. As 
we saw in section 3, the scattering amplitudes can be written as a sum of correlators in a 
grand canonical ensemble of unitary matrix models, with time setting the fugacity: 

oo „ 

F(x ,z,u c ) = J2(~ X ) N / dM u{N) \det{l-zM)\ 2i "° , (C.l) 

AT=1 

where x = 2nXe x °. Using techniques developed in J2JJ, we can show that this can be 
written as 

oo 

F(x,z,E 2 ) = J2(-x) N (mef dmHmm(e) \N) , (C.2) 

N=l 

where f(6) = — //ln(l + r 2 — 2rcos9) and \N) is the A-fermion vacuum state satisfying 

A - 1 

B n \N) = 0, \n\ > — — , 

N-l ™ 
Bt n \N) = 0, W<— 2~ • 

Here, we have expanded the fermionic fields ^(0) and ^'(9) in terms of creation and 
annihilation operators as 

V(0) = B n e in \ #+(0) = Bl n e~ in9 . (C.4) 



The expression in (|C.2 ) is still quite cumbersome because we have to sum over expectation 



values in the N fermion ground state. To obtain a more useful expression, we will add 
a factor of e~@ H and sum over all states in the Hilbert space. Then, by projecting onto 
the N fermion sector and taking the j3 — > oo limit correctly reproduces the sum over N 
fermion ground states. Thus the sum in ( |C.2| ) can be written as the f3 — > oo limit of the 
following partition function: 

Z(J3) = Tr(e-P H ef dmi W f ^ 9 \-x) N P N ) , (C.5) 
where Pn is the projection operator onto the N particle sector: 

P N = [dae ia ^- N ) . (C.6) 
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Here, N = J d9^^ counts the number of fermions. The Hamiltonian H in the iV fermion 
sector is given by 

H = H + E Q --Jde—— , (C.7) 

where we have defined Hq = | J d9 ^r^f an d E = — N ( N ^~^ _ Notice that the N 
fermion ground state \N) has zero energy. Using these definitions we can write ( |U.5| ) as 

Z(fl = f daTr(e~ 13110 el ( )(K 9 )+*<x) 9 (. ) \ e -i(a+ir)N+/3E +lnx _ ^ c _g^ 

To proceed further, we notice that the trace in this expression can be computed by 
path integral methods.: 

Tr(e- 0H O) = [ V^V^e~ s O . (C.9) 

The path integral naturally computes the expectation values of operators which are normal 
ordered, which here means that all the \l/''s are pushed to the left. Thus to compute the 
trace 

Tr(e _/3Ho e/ de * t W(/W+ iQ )*W^ (C.10) 

from a path integral, we write the operator e-^ de * ^^f^ +ia ^^ as a normal ordered 
operator: 

Then, the trace can be written as a path integral expression as follows: 

^-PHoJae^ieKm+iamo^ = |^tpvl> e / ded^lI ' t (^ + ^ +e(/(e)+^Q)5 ( ^ )) ,I ' 

= det(z^ + ^ + e lQ (l + r 2 - 2r cos 0) 

(C12) 

where the path integral is performed over periodic \I/'s satisfying \I/(0, #) = ^((3,6). This 
determinant is still non-trivial to compute and here we do not push it any further. If this 
determinant is computed, then the (3 — > oo limit of Z{(3) in (|U.8|) can give the amplitude 



Appendix D. Relation to a grand canonical ensemble of a classical gas 

It is interesting to note that the series (|5.2|) can also be rewritten as the grand canonical 
partition function of a Coulomb gas of point charges in two dimensions, confined in a circle 
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and interacting with two external point charges via repulsive Coulomb interactions [5S 
The point charges on the circle interact via a two body potential 

V(t h tj) = -ln|e^ -e**'| . (D.l) 

The vertex operators at Z\ = and = r correspond to two external point charges at 
= oo and ljz\ = 1/r of strength —a, acting on the charges on the circle by the 
potential 

V(tj) =a Hl-* e w *| =«ln|l-re^| . (D.2) 

a=l,2 

Now, we can introduce a temperature 1/(3 = 1/2, and the coefficients J n ( |5.3[ ) can be 
interpreted as the partition function of a gas of n charges interacting with the two external 
ones at temperature T = 1/2: 

Z N (T) = J N = J] J ^ e-^I- ^^) + ^.<. ^C*-'*')) . (D.3) 



If we formally introduce a fugacity z = = — 27r\e x , the infinite series ( |5.2| ) becomes the 
grand canonical partition function: 

oo 

Q(T, z) = F(r, a; x°) = £ z N Z N (T) . (DA) 

Note that for A > 0, the fugacity has a "wrong" negative sign. Further, eventually we will 
want a = iu>2 , so that the strength of the charge at Zi = 1/r becomes imaginary. 
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